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Abstract 

One-loop corrections to kink masses in a family of (l-l-l)-dimensional field theoretical 
models with two real scalar fields are computed. A generalized DHN formula applicable to 
potentials with and without reflection is obtained. It is shown how half-bound states arising 
in the spectrum of the second order fluctuation operator for one-component topological kinks 
and the vacuum play a central role in the computation of the kink Casimir energy. The issue 
of whether or not the kink degeneracy exhibited by this family of models at the classical level 
survives one-loop quantum fluctuations is addressed . 



^ Introduction 

. P.' 

BPS states arising both in extended supersymmetric gauge theories, [1], and string/M theory, [2], 
play a crucial role in the understanding of dualities between the different regimes of these systems. 
In this framework, domain walls appear as extended states in N=l SUSY gluodynamics and the 
Wess-Zumino model, [3]. It is desirable to compute the one-loop quantum corrections, AM, to 
the masses of these new entities. The huge number of fields involved in these theories, however, 
renders the task impossible. In a search for inspiration about this problem, study of quantum 
corrections to the masses of (l+l)-dimensional real solitons has been reignited in recent years, 
both in a supersymmetric and a purely bosonic framework. 

This topic was first addressed in the classic papers of Dashen, Hasslacher and Neveu, [4], 
and Faddeev and Korepin, [5] . The authors treated the purely bosonic A[</)]2 and sine-Gordon 
models for a single real scalar field. In the first model, the quantum correction to the kink mass 
was established to be AM = hm — The response given in those papers is currently 

accepted to be correct. In the eighties, the supersymmetric extension of these theories was studied 
widely in references such as [6, 7, 8, 9, 10]. The main concerns were the quantum correction, AM, 
to the SUSY kink mass, and whether or not the quantum Bogomolny bound was saturated. By 
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the end of the past century consensus about both questions was reached, mainly based on the 
papers of the Minnesota, Stony Brook- Vienna and M.I.T. groups: [11], [12] and [13]. In [11] using 
powerful supersymmetric techniques a new anomaly in the central charge was found, following a 
conjecture in the second paper of [12]. In [12] a profound analysis of the bosonic issue was achieved. 
These authors carefully distinguished between the outcomes obtained using two different kinds 
of cut-off regularization procedure. They found that only the regularization method based on a 
subtle cut-off in the number of fluctuation modes flts in smoothly with super symmetry, precisely 
the method that leads to the same result as in the computation performed by DHN/FK for the 
bosonic fluctuations. We shall refine the mode number regularization procedure in a way suitable 
to be applied also to potentials with reflection, carefully approaching the problem from the ID 
Levinson theorem as developed in Reference [14]. Dirichlet boundary conditions will be imposed 
to manage a discrete spectrum; a cut-off in the energy will be considered, and a finite number 
of modes near threshold must be subtracted in such a way that the numbers of modes counted 
around the kink and the vacuum arc identical; after this, the continuous limit is taken. 

Graham and Jaffc [13] also obtained this response, using techniques based on continuous phase 
shift methods. Moreover, they extended the analysis to the trivial sector in order to compute the 
one-loop correction to the static energy of one (far separated) kink-antikink configuration ^. Here 
the authors also notice one important point: half-bound states in the spectrum of the second 
order fluctuation operator around the static conflguration require a sharper treatment. Weak 
fluctuations around the A (0^)2 and sine-Gordon kinks are governed by Schrodinger operators which 
involve potentials without reflection. Thus, the kink zero-point energy receives a contribution from 
half-bound states that is exactly cancelled by the subtraction of identical contributions from half- 
bound states to the vacuum Casimir energy. In general, this is no longer so and extreme care is 
needed in dealing with the non-pairing of half-bound states to the kink and the vacuum sectors. 

In both [15] and [16], generalized zeta function and heat kernel regularization methods have 
been used to compute quantum corrections to the mass of SUSY (the flrst paper) and bosonic 
(the second one) one-component kinks. This latter procedure directly uses information coming 
from the potential - no need to unveil the spectrum - of the Schodinger operator and one skips 
the subtleties posed by half- bound states. In fact, the heat kernel high-temperature expansion, 
see [17], allowed us to express the one-loop correction to the A(0^)2-kink mass as an asymptotic 
series with a relative error of the 0.07 %. 

The main merit of this approach is the breaking of the number of fleld components barrier. So 
far, only kinks with a single non-null component were susceptible of being treated semi-classically; 
the difficulty lies in the study of the spectrum oi k x k non-diagonal matrices of differential 
operators. In Reference [18], however, the one-loop correction to the mass of two-component 
topological kinks in the celebrated MSTB model -[19] - has been given as an asymptotic series, 
starting from the heat kernel expansion. 

In order to be complete we also mention that there are other two interesting regularization 
methods to calculate corrections to kink masses. First, see [20], a local mode regularization - the 
cutoff reflecting the spacial variation of the kink - conceptually improves the mode number method 
and leads to the correct answer. Second, see [21], the dimensional regularization procedure is the 
appropriate method and has been successfully used when the kink is embedded in a domain wall. 

In this paper we shall discuss a model encompassing two real scalar fields proposed by Bazeia 

-^The kink-antikink system was first correctly analyzed by Schonfeld [9]. 
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and coworkers, see [22] , where the next level of complexity in computing quantum corrections to 
kink masses arises. The main novelty with respect to the MSTB model is the existence of classically 
degenerate kink famihes, see [23] and [24]. The importance of the model hes in the fact that 
Shifman and Voloshin [23] have shown this system to be the dimensional reduction of a generalized 
Wess-Zumino model with two chiral super-fields. Thus, in this framework the kink solutions turn 
into BPS domain walls of a effective supersymmetric theory. An explicit demonstration of the 
stability of some of these solutions is presented in ref. [25] using techniques of SUSY quantum 
mechanics. More recently, see [26] , it has been shown how to modify the supersymmetric version 
of this model to make the system compatible with local super symmetry. The kinks of the (1+1)- 
dimensional model become exact extended solutions of A/" = 1 (3+l)-dimensional supergravity 
with the local superpotential proposed by Eto-Sakai. Study of the effective dynamics arising from 
quantum fluctuations around these kinks is of great interest because they become exact domain 
walls (two-branes) in J\f — 1 supergravity. 

Whether or not the BPS kinks remain degenerate in mass -or the BPS domain walls in surface 
tension- after taking into account one-loop quantum corrections is the main concern of this inves- 
tigation. A clue to answering this question is offered in our paper [27]. The low-energy dynamics 
of BPS kinks is shown to be determined by geodesic motion in the kink moduh space for a special 
value of the coupling constant. Bohr-Sommerfeld quantization of this adiabatic evolution amounts 
to treating the Laplacc-Bcltrami operator of the metric inherited from the kink kinetic energy as 
the Hamiltonian. The kink moduli space is the half-plane and the metric is flat: Therefore, the 
Hamiltonian is the ordinary Laplacian for an appropriate choice of coordinates and it seems that 
the quantum effects do not lift the kink degeneracy, at least in this quantum adiabatic limit. 
We shall show, however, that these expectations are not fulflUed at the semi-classical level and 
repulsive forces between separate lumps arise. 

The organization of the paper is as follows: in Section §2 we introduce the model and briefly 
describe the structure of the conflguration space and the rich variety of kinks. Section §3 is 
devoted to computing the quantum correction, AM, to the one-component topological kink TKl 
by means of a generalized DHN procedure applicable to potentials with and without reflection. The 
contribution of the half-bound states as a function of the coupling constant is carefully analyzed. 
In this Section the one-loop corrections to the mass of the TKl kinks are also estimated, using the 
zeta function regularization method and the high-temperature expansion. Comparison between 
the approximate and exact results for AM(TKl) serves as an evaluation of the error. Starting 
from numerical solutions. Section §4 offers a computation of the semi-classical masses for the 
whole BPS kink variety at three special values of the coupling constant. If cr = 2.5 and a — 1.5 we 
flnd a rapid convergence of the asymptotic series, although violation of the classical degeneracy 
is suggested by our results. The anomalous kink degeneracy is analyzed in connection with the 
existence of BPS link kinks. For a = 2, the system becomes equivalent to two decoupled A[0]2 
models and the numerical method provides results in agreement with our previous work in [16]: 
in this case the kink degeneracy still holds at the one-loop level. We also compute in this Section 
the one-loop corrections to the mass of the BPS link kinks for the value o" = 2 of the coupling- 
constant in order to easily show that the kink sum rules hold at the semi-classical level. In Section 
§5 the dependence of the one-loop kink mass formula on the parameter c is studied from a general 
point of view, relying once again on asymptotic methods (the high-temperature expansion) . The 
model also has non-BPS kinks. Section §6 attempts to elucidate why non-BPS kinks arise in some 
regimes of the coupling constant. In particular, study of stability via index theorems shows that 
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families of this kind of kinks exist for some critical values of the coupling constant. Finally, an 
Appendix is offered explaining the generalized DHN formula, derived using the cut-off in the mode 
number regularization procedure but also valid for potentials with reflection. 



2 The model: classical kink degeneracy 

We shall focus on the (l+l)-dimensional two-component scalar field model introduced in [22]. 
This system arose as the bosonic sector of a supersymmetric theory with superpotential: 



W(x) 



1 2 



Xi{y^)^i+X2{y^)(^2 is a two-component real- scalar field; y^, = 0, 1 , are coordinates 



^'^ space-time, and ea , a = 1, 2 , form an orthonormal basis in the the internal space : 
Wc also choose the metric g^^ = diag(l, —1) in R^'^ and a system of units where c = 1 



Here 
in the I 

4 ■ 66 = Sab 

but keep h explicit. The coupling constants thus have the following units: [A] = — M~iL~^ 
and [a^] = ML. The dimension of the scalar field is: [x] = M^L^ . 
The dynamics is governed by the action: 



Six] 



Id^X-d^X-l^W-VW 



(1) 



Introducing dimensionless fields, variables and parameters x — 2a0, — fx^ — ^ a' 
we obtain expressions which are simpler to handle. 

For static configurations and the dimensionless superpotential — 2 — -|- f^i^l) 
the energy functional reads: 



%i,X2]=4a=^A£:[0i,02] 



dx 



(2) 



It is a function of the unique classically relevant coupling constant a and also depends on the 
dimensionless potential term ?7(0i, ^2) = | (40i + 2cr0| — 1)^ -|- 2(T^0^0|. 

There are four classical vacuum configurations in this system, organized in two orbits of the 
internal parity symmetry group G/ = Z2 x Z2 generated by the transformations tti : (0i,02) 
(—01,02) and 772 : (0i,02) (0i, —02)- The vacuum moduli space, however, includes two points 
if CT > 0: 

Jvi = — = pomt U pomt 
Gi 



Gj_ 

r(i) 



U 



(pv,±{x,t) 



(2) 



1 ^ 

= 62 



'2a 



Here, i/j"'' is the sub-group of Gj that leaves the (py± point invariant. If a < 0, the (py± is lost 
and the vacuum moduli space has only one point in this regime. 
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Figure 1: The potential —U{(j)i, 4)2) for ct < (left) and for u > (right) 



We shall restrict ourselves here to the a > regime, because a richer manifold of kinks arises 
in this range. The space of finite energy configurations C is the union of sixteen disconnected 
topological sectors in this case. If a, 6 = 1,2 we have that: 



C — Ua^b ^fi^b) '^a^b ^fa^b) 



with 



In the Cf^^ non-topological sectors a symmetry breaking scenario develops: the original Gj sym- 
metry of the action (1) is broken to the Hf"^ sub-group by the choice of |0y±) as the "zero-particle" 
state in the Cf^-, sector. 

(a, a) 

The small fluctuation solutions f]a(t,x) = 0^ + Sri(t,x) of the field equations 



-VU 



around the vacuum points and V.^ can be respectively expanded in terms of the eigenfunctions 
of the Hessian operators: 

Therefore, the one-particle states have respectively (dimensionless) masses: m\{y^) = 4, m^Vf^) = 
and mliV^) = mliV^) = 2a. 

The most important feature of this system is that there exists a degenerate in energy family 
of BPS kinks in the C^^-^ topological sectors for any value of a > 0. BPS kinks are finite energy 
solutions of the first-order system of differential equations: 

~ ~ ± 20^ 0-02 - o -r- = = ±20-0102 ■ (3) 



dx 901 \ 2 J dx 

In both References [23] and [24] it has been shown that the BPS kinks belonging to the C"^^^^ 
topological sectors are in one-to-one correspondence with the curves in the internal space: 
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2(1 -cr)^^ 4 2(7 

where c e (— oo, C^) is an integration constant and there is a critical value = \-^{2a 
a = l, 



^+ log 1021) = ^ 



c e (-oo,c^) , = -1 + ln2. 



One sees in Figure 1 (right) that there are two maxima of —U{(j)i, ^2) with the same height. Kink 
solutions which go from one maximum to the other depend on a parameter c which measures 
whether the particle moves through the bottom of the valley, or more along the sides on the curve 
(4). There is a critical value of c where the particle moves as high as possible; increasing c beyond 
this critical value the particle crosses the mountain and falls off to the other side. 

All the BPS kinks belonging to these families have the same (dimensionless) energy: £(BPS) — 
\W{(l)y±) — W{(l)y^)\ = |. We shall denote these kinks as TK2(c) because from (4) one sees that 
the two components of the field are non-zero for generic values of c; also, all the kink orbits in (4) 
start and end in vacua belonging to the same Gj-orbit, and are accordingly classified as "loop" 
kinks. 

Plugging (4) into (3) one reduces the solution of the (3) ODE system to a single quadrature. 
The explicit integration can be performed analytically only for certain special values of a, see 
References [23], [24]. 

Exactly for c = —00, one-component topological kinks - therefore termed as TKl- arise. If 
c = —00, for any a G [0, 00), we obtain the curve 02 = in (4) and find the BPS kinks in the 
^(1^) sectors of configuration space: 

1 ^ 
<pTKi{x) = (— l)";^ tanh(a; -|- a)ei a e Z/Z2 

a e R is another integration constant that fixes the center of the kink. 



3 One-loop correction to the mass of one-component topo- 
logical kinks 

In this Section we shall compute the one- loop quantum correction to the classical mass of TKl 
kinks. Because the matrix differential operator governing second-order fiuctuations around these 
kinks is diagonal, 

/ -^+4-6sech2x \ . . 

I ' -^ + <72-a(<7 + l)sech^x; 

this task can be performed using both the DHN formula and the asymptotic method. Note also 
that the one-loop correction to the TKl kink mass is the sum of two contributions: (1) AM(/Cii) 
is the one-loop correction due to the the tangent fiuctuations to the TKl kink governed by /Cn. 
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(2) AM{IC22) collects the contributions to the quantum corrections coming from the orthogonal 
fluctuations to the TKl ruled by /C22- Therefore, we write 



AM(TKl) - AM{]Cu) + AM{IC22) 



3.1 One- loop correction to the TKl kink mass: DHN formula 

In (5), both JCu and /C22 are particular cases of Schrodinger operators of Posch-Teller type 



+ Co- 



Uo 



cosh X 



whose eigenvalues and eigenf unctions are known [28] . Since /C22 is a Schrodinger operator involving 
potential terms with reflection {a ^ N) and without reflection (a e N), we use the generahzed 
DHN formula for the /l-correction to the mass of one-component topological kinks: 



AM{]Caa) = 



hm 



i=0 



Va ^ 1 d5a{q) r^— ^ {Vgajx)) 

2^ 



+ 



+ hm 



{Vaaix)) 



dk 



(6) 



where 



Vaa{.x) 



[X] 



and (■) stands for the expectation value: {Vaa{,x)) = dxVa, 

Formula (6) has been derived in the Appendix collecting previous work in this topic. We set 
as m = aA the parameter with dimension of that will be used to fix the real dimension of 
each observable in the system. Note that as a consequence of our choice of dimensionless variables 
there is a global factor of \/2 in formula (6) with respect to the analogous formula used in [16]. 
We now pause to explain the origin of the different terms: 

1. The first line accounts for the zero-point energy of the quantum kink measured with respect 
to the zero-point energy of the vacuum. Eigenstates from both the discrete and continuous 
spectrum of ICaa contribute. The highest bound state contributes with a weight s/, which 
is 1 if uji does not coincide with the bottom of the continuous spectrum (threshold). If ui 
is buried at threshold, the corresponding eigenstate is a half-bound state and = |. To 
deduce formula (6) the density of states in the continuous spectrum is given in terms of the 
phase shifts Sa{q) of the scattering waves through the potential: — Vaa{x). A cut-off in 
the number of modes (sec Appendix for a complete derivation of formula (6) including the 
case of potentials with non-zero reflection coeflicicnts) has been used for renormalizing the 
zero point energy. Note also that the subtraction of ^ amounts to taking into account the 
contribution of the half-bound state in the vacuum to the vacuum zero point energy (see 
[13] and the Appendix). 

2. The second line takes into account the contribution of the mass renormalization counter-term 
to the Casimir energies of both the kink and the vacuum. 
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3.1.1 Tangent fluctuations: AM(/Cii) 

Spcc(/Cii) inchidcs two bound states, one half-bound state (in the above notation 1 — 2), and 
scattering eigenfunctions characterized by the phase shifts: 

3q 

— — 2arctan; 



2-g2 
Therefore: 

Spec(/Cn) = Wo = 0} U {uJi = 3} U {uj2 = 4},^^i U {q" + 4},eK 
Substituting all this information into formula (6) we obtain: 



AA^/^TKu VShm hm 3y^^niq^ + 2) , 3hm dk Hm f 



dx6 sech^ x 



= hm 



3 



2^3 TT 



Except for a global factor -\/2, the contribution of the tangent fluctuations is the same as the 
contribution of all the fluctuations to the kink in the 0^ model; note that there are two bound 
states with eigenvalues of and 3, giving tangent fluctuations to the TKl kink identical to those 
arising in the A(0^)2 model. The global factor is due to the different choice of the parameters in 
the potential. 

3.1.2 Orthogonal fluctuations: AM(/C22) 

The scattering in the potential well - Cq = cx^, uq = a{a + 1) - of /C22 is not "refiectionless" if a is 
not an integer. There are in general even and odd phase shifts 



'™ 4 \Re(T(q)±R{q})J 



to be read from the transmission and reflection coefficients 

^ r{a + l-iq)T{-a -iq) ^ r(a+l- iq)T{-a - iq)T{iq) 

r{i - iq)r{-iq) ' r(i + <7)r(-<7)r(-ig) 

RecaUing that 

e^'^t ^ T{q) ± R{q) (7) 
T{q) and R{q) are obtained from the asymptotic behaviour of the scattering eigenstates: 

iP,ix) = Ne"^\F,[-a, a + l,l-zq, , 

where 2-^1(0, b, c; d\ is the Gauss hypergeometric function . 

If cr = Z G N is a natural number, R{q) = - r(— cr) has a pole in this case-, 62 {q) — (^2^(?), 
and the total phase shift ^2(5) = S2{q) + ^2 (?) i^- 



S2iq) ^ iarctan ( l^nU^^i^^ - H - + ^ - n)) ^ 
2 I Rento(?^ - a - ny + 2iq{l - n)) ' 
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The number of bound and half-bound states of Spec(/C22) is equal to /[cr], I[a] being the integer 
part of a. The corresponding eigenvalues are cui = i{2a — i), see [29, 28], whereas the bound state 
eigenfunctions read: 

N 

V'i(a^) = ■^^^^j^^y(— y2-^i[-«,2a - i + 1,C7 - i + 1, 1(1 + tanhx)] , i = 0, 1, • • • /[cr] 

The half-bound state arises only in the case a = I[a] = I eN. Some explanation about this issue 
should be provided: The "first" eigenfunction in the continuous spectrum 

ipg=o{x) = N2Fi[-a, a + 1,1; ^(1 + tanhx)] 

does not belong to the Hilbert space because limx^±^ipq=o{x) = c±x, except if uq = l{l + 1). In 
this latter case, -01=; = ipq=o , 

lim 2Fi[-l, / + 1, 1; i(l + tanha;)] = c± 

x-^±oo 

and hence this is the half-bound state in the spectrum of 1C22- 
In sum, we have that: 



Spec (/C22) 



Ui=o,i,...,/HK = ^(2(7-«)}U{g2 + a2},gR+ if cr^N 

Ui=o,i,...,a-iK = «(2(T-i)}U{a;,=^ = (72}^^^^ 1 U{g2 + (T2}geu+ if aeN 



Therefore, in the subspacc of C^\-j orthogonal to the TKl configuration more bound states, de- 
pending in number on the value of cr^, appear as eigen-fiuctuations of the Hessian operator. The 
corresponding eigenvalues are semi-definite positive for any value of sigma. The TKl kink is thus 
stable independently of the parameter. An important point is that there exists a second zero 
mode bccaTisc the first eigenvalue of /C22 is zero for all a^. This zero mode obeys neutral equi- 
librium fluctuations in a continuous family of kinks with the same energy as the TKl kink and 
second non-zero component, see [24]. 

Formula (6) demands the computation of the derivative of the phase shift with respect to the 
momentum q, using (7) can be performed explictly: 



dq ~ 2 



,2stdf±^^-m^de^ 
dq dq 



= 2Remq)-il;{-a + iq)] + ——^ — — . (9) 

2 smh 7rg esc iTxa -\- tan ttct 



Plugging all these expressions into (6) we have that: 



AM(/C22) 



hm 

~Y 

hm 



o-(l+cr)\ a{a + l) 



a-l 



" ^ poo 

^^/^{2a-i)+- J dq 



,i=0 



d52{q) 
dq 



cr(l + 0-) \ cr(cr + 1 



TT 



TT 



if (T^N 

if ctGN .(10) 
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We are now interested in describing separately the behaviour of the contributions from the discrete 
spectrum of the operator /C22 and from the integral in (10) to AM(/C22)- We define 

AM(4?) = Ev^(^ , = f^y?T^+ 4=) (11) 

to emphasize the different contributions. 

In AM()C^), the contribution of the discrete spectrum of K.22 to the one- loop correction of the 
TKl kink mass, wc notice the following. First, the greater the value of a, the greater the number 
of bound states in Spec(/C22) is. There is a jump in AM(/C5J) when a crosses an integer value, 
sec Figure 2. Second, a half-bound state in Spcc(/C22) arises exactly at the values a = I E N, 
which contributes to AM(IC'^) with a | weight with respect to the contribution of any other 
state. One would think from these considerations that AM(/C22) is a non-continuous function of 
the coupling constant a because of the jumps in AM{]C^). On the other hand, we know that 
the TKl solution depends smoothly on the value of cr, such that a smooth response of AM(/C22) 
would be expected. 



15 
12.5 
10 
7.5 
5 

2.5 




Figure 2: Contribution to the one-loop correction from the discrete spectrum 0//C22 (left) Behaviour of 
the integrand f{q, a) in (10) (right) 



The analysis of the second contribution is also subtle. The main points are as follows: First, the 
descriptive behaviour of the integrand /(g, cr) depends only on the fractional part of o", following 
the pattern shown in Figure 2 for the range a e (1, 2). Second, asymptotically we have that: 

lim /(g, cr) = , V(T . 

q— >oo 

Moreover, 

hm /(g, a) = 1 + a(l - 2^e - V'(-fT)V^(l + a)) , 

where 7^ is the Euler Gamma constant, and ^^{z) is the Digamma function. Note that /(g, a — I) 
has a pole at g = 0. Because fiq.cr) is wcll-bchaved for any other value of g if cr is non-integer, 
the integral in (6) is well defined and a finite answer for AM(/C22) is obtained, see Figure 2. The 
integral becomes improper if cr = / e N: in this case it is only defined by taking the limit g — >^ 
appropriately. A miracle happens: the integral of the function /(g, cr) produces a jump in this 
contribution when cr crosses an integer that exactly cancels the jump in AM{1C^^) due to the 
appearance of a new bound or half-bound state!, AM(/C22) is a continuous function of cr. 



10 



Nevertheless, the integral cannot be evaluated analytically except for integer values of a. We 
use numerical analysis to compute AM(/C22) by means of formula (6). The numerical computation 
confirms that AM(/C22) is smooth in cr. The contributions from the discrete and continuous 
spectrum balance each other to obtain a smooth behaviour in a. We stress that this balance is 
made possible by a proper counting of the half-bound states in both the vacuum and TKl kink 
Schrodinger operators. 

In the next Table we offer the data on the sum of the tangent and orthogonal contributions to 
the TKl quantum mass correction, AM(TKl), for some values in the range a e (0.4, 3.3), which 
are also depicted in Figure 3. 



a 


AAf(TKl)/ftTO 




a 


AM(TKl)/ftm 


0.4 


-0.799335 




1.4 


-1.11618 


0.5 


-0.829892 




1.5 


-1.15057 


0.6 


-0.860369 




1.6 


-1.18559 


0.7 


-0.890955 




1.7 


-1.22128 


0.8 


-0.921788 




1.8 


-1.25765 


0.9 


-0.952966 




1.9 


-1.2947 


0.99 


-0.981384 




1.99 


-1.32865 


1.00 


-0.984565 




2.0 


-1.33251 


1.01 


-0.98775 




2.01 


-1.33627 


1.1 


-1.01664 




2.1 


-1.37094 


1.2 


-1.04925 




2.2 


-1.41013 


1.3 


-1.08242 




2.3 


-1.45005 



(T 


AM(TKl)/ftm 


2.4 


-1.4907 


2.5 


-1.53212 


2.6 


-1.57427 


2.7 


-1.61717 


2.8 


-1.65316 


2.9 


-1.70527 


2.99 


-1.74592 


3.0 


-1.72309 


3.01 


-1.75503 


3.1 


-1.79644 


3.2 


-1.84319 


3.3 


-1.89071 



3.2 Asymptotic series for the TKl semi-classical masses 

The evaluation of one-loop corrections to the masses of other (two-component) kinks which are 
classically degenerate in energy with the TKl kink is another task to be addressed in this work. 
For this purpose the DHN formula is of no use, because our knowledge about the spectrum of non- 
diagonal matrix Schodinger operators is grossly insufficient. Alternatively, we can use a formula, 
derived in References [16, 18], that expresses the one- loop correction to /c-component kink masses 
as an asymptotic series. If we denote as /C the Hessian operator around a given loop kink K, the 
formula for the one-loop correction to the masses of topological kinks derived in Reference [18] 
reads: 



3 



AM(/C) = /im[Ao + D, 



D. 



no 



2^ 

2 no-1 

-EE 

a=l n=2 



8n 



,2n-2 



UqEN 



(12) 



A lightning summary of the content of formula (12) is as follows: j = dim Ker(/C) is the number of 
zero modes in the spectrum of /C. [a„]aa(^) ai'e the coefficients of the high- temperature expansion 
of the heat function associated with the heat equation: 



ab 



7[n — 1,^^] are incomplete Gamma functions, see [30], and Va are the masses of the fundamental 
quanta at the point of the vacuum moduli space determined by the loop kink K. The divergent 



11 



terms 



E 



Naa(^)7[-l,^^a] 



E 



[ai]aa(^) 



a=l " a=l 

lacking in (12) are respectively cancelled by subtracting the zero point vacuum energy, and taking 
into account the mass renormalization counterterm, see [16] and [18]. 

Although we have computed ^M{TK1) by means of the DHN formula (6) in the previous 
Section, we now apply (12) to the TKl kink in order to test how good the approximation provided 
by the asymptotic series is with respect to the exact DHN result. 

The spectrum of /C comprises two zero modes; henceforth, j = 2. Moreover, vf = A and 
V2 — in this case, whereas the Seeley coefficients [a„(/C)]ii and [an(^)] 22 are pointed out in 
References [16, 18]. The calculations of the partial sums in (12) for several values of a and with 
Uq — 10 provide the results shown in the next Table: 



cr 


AM(TKl)/fim 




a 


AM(TKl)/;iTO 


0.5 


-0.962386 




1.1 


-1.05073 


0.6 


-0.970537 




1.2 


-1.07468 


0.7 


-0.981183 




1.3 


-1.10097 


0.8 


-0.994487 




1.4 


-1.12939 


0.9 


-1.01053 




1.5 


-1.15971 


1.0 


-1.0293 




1.6 


-1.19174 



a 


AM(TKl)/fim 


1.7 


-1.22526 


1.8 


-1.2599 


1.9 


-1.29571 


2.0 


-1.33324 


2.1 


-1.37074 


2.2 


-1.41007 



In Figure 3 the results are plotted as (white) squares and are compared with the data obtained 
from the DHN procedure, (black) dots. We can check the concordance between the two methods; 
the relative error for the answer from the asymptotic method is, for instance, 0.79 % for o" = 1.5, 
0.055 % for cr = 2.0 and 0.004 % for cr = 2.2. The greater cr, the more exact the response reached 
from (12). On the other hand, for the range cr € (0, 1) the answer is less precise (for a = 0.9 the 
relative error is 6.0 %). The reason for this becomes clear taking into account that the greater the 
value of the masses, the smaller the value of ^^"^^i'^"^ , and the more rapidly the series expansion 
in formula (12) converges. 
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Figure 3: One-loop correction to the one- component topological kink (TKl) mass in units of hm. •, 
DHN formula . asymptotic series. 

It is interesting to note that the answer obtained by means of the asymptotic series is a 
continuous function of a. This method, based in the heat kernel expansion and the associated 
generalized zeta function, leads directly to the correct result, without needing to carefully balance 
the contributions from the discrete and the continuous parts of the spectrum. The reason is that 
the Seeley coefficients are defined in terms of the potential (and its derivatives), which encodes all 
the properties of the spectrum. 
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4 Semi-classical masses of kink families 



So far we have only computed the one- loop correction to the classical mass for the TKl kink- 
using either the (exact) DHN formula (6) or the (approximate) asymptotic expansion (12). The 
model, however, exhibits a continuous family of kinks which have the same classical energy. Our 
goal in this Section is to compute the one-loop correction to the kink mass for each member of 
the family and to analyze the fate of the kink degeneracy in the semi-classical level, using formula 
(12) with the help of some numerical analysis. 



4.1 (7 = 2 

In this case, it is possible to find explicit analytical expressions for the solutions corresponding to 
the kink orbits (4), see [24]: 

r , ,T (-1)° sinh2(a; + a) ^ (-1)^ ^ , , 

(^TK2[x; a, b] = ^ \ ^ ei + ^ ^ 62 , 13 

2 cosh 2{x + a) + 2 cosh 2{x + a) + 

where a e R, 6 = ~ ^c'^-ie ^ (1' °°)' '^iP ^ Z/Z2. The family of Schrodinger operators ruling 
the small fluctuations in the original fields is in this case 



_ d? _l_ f^ smh^ {2x)+b^ -1 _ r, r) /jo 1 sinh(2a:) 

(cosh(2a:)+6)2 ^ U I (cosh(2x)+fe)2 

10,/;,2 1 smh(2a:) , n Sin}o?{2x)+h'^-l 

IZV O 1 (cosh(2i;)+6)2 dx2 "T U (cosh(2a;)+6)2 



The asymptotic method can be apphed by substituting the corresponding expressions for this 
case in formula (12): j = 2, t;^ = 4, a = 1,2, and the potential terms of the above operator are 
required in the definitions of the Seeley coefficients [a„]aa(^)- This calculation has been carried 
out for some kinks of the family (13) with integration constant c in the range c G [— 30,c'^), see 
Table and Figure 4. In sum, the approximate value of the one- loop quantum mass correction to 
all these solutions is — 1.33280 ± 0.00001 in units of hm. Thus, the statement that the degeneracy 
of the classical mass to the kink family in the case cr = 2 is preserved in the quantum framework 
is highly accurate. 

For this value of a we can easily prove that the statement is completely exact. A rotation of 
45° in R^, ei = -^{^i + ^2), &2 = 75(^1 ~ ^2), shows that the system is non-coupled. Writing 

(j) = ipiSi + ip2£2, we have that: 

We can write the degenerate kink family in the alternative form: 
_ (~1)" 

(/>TK2*[x;ai,a2] = ^ tanh(a; + ai)£i H ^ tanh(a; + 02)^2 , 

2v2 2v2 

in terms of the new parameters ai, 02 G {—00, 00). Note that the TKl kinks correspond to Oi = 02. 
In these variables the Hessian operator is diagonal for any member of the TK2 family: 



/C(ai, 02) 



, 

dx^ cosh^(a;+ai) 

(f^ I /I 6 

dx^ ~ cosi?(x+a2) 
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Therefore, both the spectrum of /C(ai,a2) and the one-loop correction to the kink masses are 
independent of ai,a2 G (-00,00): AM(TK2*[ai, aa]) = hm{^ - f ) ~ -1-33281. The kink 
degeneracy is not broken by quantum fluctuations at the one-loop level; this result would be 
expected from general arguments because the family of second-order fluctuation operators above 
is isospectral. 

4.2 Numerical approximation 

The problem that we face is the lack of explicit analytical expressions of the kink solutions for 
generic values of cr; only the kink orbits (4) are available. However we can solve the first-order 
equations (3) by standard numerical methods. We set the "initial" conditions: 



The rationale behind this choice is as follows: (1) for any kink solution, (j)i{x) has always a zero. 
Translational invariance allows us to set the zero at x = 0. (2) To ensure that we will find a 
numerical kink solution, we fix 02(0) in such a way that (4) is satisfied for a given value of a and 
arbitrary choices of c. 

In this way the numerical method provides us with a succession of points of the kink solution 
generated by a interpolation polynomial. The kink polynomial is then used to calculate the 

[an]aa(A^) cocfficicnts of formula (12) approximately and hence the one-loop correction to the 
kink mass. A subtle point needs to be clarified: the coefficients are defined in terms of the field 
solutions and their derivatives. Whereas it is all right to describe the kink solutions approximately 
by interpolation polynomials, taking "derivatives" of such discrete configurations of points would 
induce important errors. Fortunately, we can use the first-order equations (3) to express all the 
field derivatives as polynomials in the fields. 

This procedure was carried out for values of a in the (1.3,3.3) range with Aa = 0.1 finding 
in all cases a similar pattern in the behaviour of the AM{K{c)). In the Figure 4 we show the 
results obtained using this method for several values of c in the range c G (— 30, c^) in the cases 
(7 = 1.5 and a — 2.5. A first observation is the perfect agreement reached between the numerical 
approximation and the exact result of the a = 2 case. Next, we notice that the greater the value 
of cr, the better the convergence of the asymptotic expansion to the exact value for kink solutions 
with c sufficiently distant from . We check that the kink degeneracy is always maintained in a 
wide range of c, but starts to fail when c approaches c^. 
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fT=1.5 (7 = 2.0 (7 = 2.5 



c 


AM 


c 


AM 


c 


AM 


-30 


-1.16009 


-30 


-1.33281 


-30 


-1.52784 


-27.5 


-1.16017 


-27.5 


-1.33281 


-27.5 


-1.52782 


-25 


-1.16128 


-25 


-1.33281 


-25 


-1.52780 


-22.5 


-1.16042 


-22.5 


-1.33281 


-22.5 


-1.52778 


-20 


-1.16061 


-20 


-1.33281 


-20 


-1.52774 


-17.5 


-1.16088 


-17.5 


-1.33281 


-17.5 


-1.52769 


-15 


-1.16128 


-15 


-1.33281 


-15 


-1.52760 


-12.5 


-1.16193 


-12.5 


-1.33281 


-12.5 


-1.52744 


-10 


-1.16313 


-10 


-1.33281 


-10 


-1.52711 


-7.5 


-1.16597 


-7.5 


-1.33281 


-7.5 


-1.52626 


-5 


-1.18205 


-5 


-1.33280 


-5 


-1.52285 


-4.6801886 


-1.24345 


-4.001 


-1.33280 


-4 


-1.52168 


4.68018860186678332 


-1.25103 


-4.00001 


-1.33280 


-3.97 


-1.52915 



-3.96594571 -1.55402 
-3.96594570565808127 -1.56127 
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Figure 4: The One-loop Quantum Mass Correction in the cases a = 1.5, a = 2.0 and a = 2.5. 



Thus, the previous Table and Figure suggest the survival of kink degeneracy at the one-loop 
level up to values of c close to . A more precise idea about what is going on in the vicinity of 
is given by the following observations: 

1) We first notice that the c = kink orbits belong to the ^^"^^^(0 7^ b) topological sectors 
for any value of a. Thus, this kind of kink is called a "link" because the associated orbits link 
different points in the vacuum moduli space. The trajectories of these link kinks enclose the orbits 
of the kink family studied previously in Section §4. For a = 2, these TK2L topological kinks of 
link type are: 

0tk2l(x) = -^^^(1 - (-1)^ tanh(a; + a))ei + ^^^(1 + (-1)^ tanh(x + a))e2 
with a G M, a, /3, 7 G Z/Z2. The link topological kinks for cr = | are 

0tk2l(x) = (1 - (-l)'^tanh^) ei + ^^^1 + (-1^ tanh 

with a,/3,7 G Z. There are no analytical expressions available for other values of a. 
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2) Second, we realize that because the existence of these enveloping kinks the energy density 
of any kink obtained for c near is formed by two lumps. We explain the general situation in 
the cr = I case. Even though the asymptotic method works poorly for this value of a (with errors 
higher than a 18%), the analytical information allows a full description of the situation. The TK2 
kink family is given by: 



1 r w (~^)° sinh(x + a) ^ / h 

(pTK2[x] b] = — ei + (-Ij'^W — ^ ^ I 

/ cosh(x + a) + y cosh(x + a) + o 

where a G M, 6 = :^7f^ G (0, oo), and Z/Z2. 

In Reference [27], we showed that the splitting into two lumps starts at c = or 6 = 1. In the 
previous subsection we have observed through numerical computations for any value of cr 7^ 2 that 
the departure from the TKl quantum correction also starts at the value of c where the splitting 
begins. Use of the c parameter is necessary to implement the numerical method, but unsuitable 
for discussing this phenomenon. Note that c G (0, c*) is tantamount to 6 G (l,+oo). & is a 
measure of the distance between the two lumps, although a highly non-linear one, see [27]. Thus, 
the breaking of the degeneracy is noticeable when the two lumps are further apart. The induction 
of this repulsive force by quantum effects can be better understood by looking at the family of 
the TK2 Hessian operators: 



d? I 3fe I 6 sinh^ x 2 6\/6sinh2; 

dx^ coshx+fe (cosha;+b)2 (cosh x+h) i 

6\/fcsinhx d? _i_ 3 b i_ 3 sinh^ x 1 



(coshx+f))T "^^^ 2cosha;+6 4(coshx+f))2 2 

In the Figure 5 we show a plot of the diagonal components of the potential for several values of c. 




Figure 5: Diagonal components of the potential for c=-30, c=-l, c=0.1, c=0.245 and c=0.249. 

There is no modification of the potential between c = —30 and c = — 1 (the off-diagonal com- 
ponents are odd functions and the area enclosed by them is zero); there must be kink degeneracy 
in this range of c. Starting at c = 0, the second component turns from potential well to barrier. 
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inducing a (still) very weak repulsion. The closer the value of c to = | , the stronger the re- 
pulsion coming from the second component, whereas the first component develops a double well, 
clearly arising from the meiosis into two lumps. Therefore, the spectrum of /C(6) is completely 
different in the regimes b € (0, 1) and b £ (10^, oo), with a transition region in between. We stress 
again that the picture is similar to this for any a ^ 2: the classical degeneracy starts to fail for the 
value of c where two peaks in the energy density appear: quantum fluctuations induce repulsive 
forces between the two basic lumps, which are constituents of a TK2 kink. The data in the Table 
also show that the rate of change of the quantum correction diminishes when the two lumps are 
extremely far apart, in a range of c where there are no significant changes in the structure of /C(c). 

We stress that the breaking of the kink degeneracy due to quantum effects suggests the pos- 
sibility of a similar phenomenon between higher dimensional topological defects as vortices and 
monopoles; it is tantahzing to think that the method here developed could also be apphed to study 
the one-loop effects on the moduli spaces of vortices and monopoles. Nevertheless, one can argue 
the survival of kink degeneracy to quantum fluctuations in the fully super symmetric system. The 
reason is the saturation of the BPS bound at the quantum level in A/" = 1 Wess-Zumino d = 2 
supersymmetric models due to equal anomalies in the energy and the central charge; see Reference 
[31] for a recent and elegant proof in systems with only one chiral superfield. In our case, two 
chiral superfields, the anomaly is given - up to some constant- in terms of the Laplacian of the 
superpotential: 

|Aiy(0(oo)) - AW^(0(-oo))| = 4 + 2(7 

a quantity independent of c. Therefore, the c-dependence in the energy coming from bosonic 
fluctuations must be exactly cancelled by taking into account the fermionic fluctuations with 
SUSY preserving boundary conditions. 

We end this Section by exploring some points about the one-loop correction to the TK2L link 
kinks. If (7 = 2, the TK2L Hessian/Schrodinger operator for a = is 



<f , A 3 



/C^((j = 2) = ( '^^'^ \ 



whereas ii a — ^, also for a = 0, we have: 



cosh 

-^ + 4 




_l._/^-^ + i-|ta^h( f)-i^ ^(l-tanh(f))Vl + tanh(f) 
^'^ 2^ I i^(l-tanh(f))Vl + tanh(f) + | + | tanh(f ) - 

In this last case, we would need to use the background renormalization method as developed in 
[16] to compute the one-loop correction, but the task is so involved that we only discuss the a — 2 
case. 

It is clear that SpecJC^{a = 2) = Specie, where: 

^ 

aa; cosh x 



One immediately concludes that 



AM{TK2L) = ^AM{TK1) 
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i.e. the famous kink sum rule also holds at the one-loop level. The results obtained via the 
asymptotic expansion agree with this observation. The very accurately calculated degeneracy in 
the case (7 = 2, obtained from numerical plus asymptotic computations even in the vicinity of the 
limiting value — —4 (see Figure 4), implies that the combination of two identical link kinks 
(reached when the parameter c goes to c^) has the same mass quantum correction as the rest of 
the kinks. 

As we have seen the situation is different (and to some extent surprising) for other values of 
a. The singular behaviour of the system for cr = 2 is related to the fact that the metric ruhng 
the kink adiabatic motion is Euclidean if o" = 2, see [27]; consequently, the two lumps move freely 
with respect to each other. 



5 Classical versus semi-classical kink degeneracy 

In this Section we shall address the subtle question of what happens to the kink degeneracy when 
quantum effects are taken into account from the analytical - more than the computational- side. 
Let us denote by 

(f)^^''\x; c) = 0i(x; c)ei + (j)2{x; c)e2 

the solution of the first-order equations (3) given by the orbit (4). Let us define /C(c) = Vi — V{c) 
as: 

jCfc) = + 2402[a;; c] + 4a(a + l)^l[x; c] - 2 8a((7 + 1)0102^; c] \ 

^' V 8a(a + l)0i02[:r;c] -£ + 4a{a + l)^l[x;c]+6am[x;c]-a J 

We have thus a family of Schrodinger operators governing the small fluctuations around any kink 
in the degenerate family. 

The one- loop correction to the mass of the K{c) kink is given by the formula (5) of Reference 
[18]: 

AM{K{c)) = lim [Ai£^(^)(s)-|-A2£^(')(s)] 

A,s^«^\s) = }irn^^^'^^'^±^j2Civ.Us + l) r dxVUx-,c) . (14) 

^ ^ a=l 2" 

We recall that Ai£^*^'^^(— |) is the kink Casimir energy measured with respect to the vacuum 
Casimir energy. A2£^*^^-'(— |), however, accounts for the contribution to the kink energy of the 
mass renormalization counterterm. 

The generalized zeta functions of /C(c) (with the zero eigenvalues excluded) and Vi can be 
written in terms of the heat functions via Mellin transforms: 



r(s) Jq r(s) Jo 

The derivative of Ai£^("^)(s) with respect to c is obtained immediately. 



^ -pPIC(c) 

dc 
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whereas we also have: 



dc 



V{s) 



mL „ 

- , dVaa 

ax — — 

mL OC 



If the sum of these two quantities adds to zero at the -s — > — | hmit we can be sure that the kink 
degeneracy also holds up to this order in the /i-expansion. Unfortunately, an exact computation is 
beyond our analytical capacities but we can work an asymptotic expansion of this formula along 
similar lines to those used in Reference [18] and [32]. Prom the high-temperature expansion of the 
heat kernel, 



n=0 



and bearing in mind that 
'dV 



Tr 



dc 



— oo 

oo 



ac 



X) 



dx 



dx'ti (x\^\x') (x'le-^^'^(^) \x) 
dc 



2 

a,b=l''~ 



dx — Kfc(c)ba{x, x; (3) 



we find 



OC 



oo 2 



^yy 



n=0 a,b=l 

oo 2 



dx 



dVabix 



dc 



47r ~: ^ 



n=0 a=l 



■[an]ba{x,x)e-^^^)f3^ 



where 



Therefore, 



[Cn]a = E / 

L 1 — 



~ dVab{x;c) ^ 

dx [an\ab{x, X) 

6=1 " ^ 



d_ 
Wc 



1 oo 2 „i 

V47rr[s]^^ Jo 



dBicjc) js) 
dc 



where S;c(c)(s) is a negligible error arising from the truncation of the (3 integration in (3 — 1. This 
gives: 



8 



hm 



oo 2 



n=0 a=l 



[Cn\a7[s + n+ l,vl] 
— 2(s+n+i) 



+ error 
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Also, we have that 



. r rnL l[s + ^,vl] 



and, because [ao]ab{y,y) = Kb, 



6=1 



5c 



Thus, 

dc 2 ^V47r^^(s + |) 

and finally we obtain: 



/I— 2( 
n=l a=l v47rVa 



It is a remarkable fact that the pole arising in lim^^_i Ai£^^"^)(s) appears in the n = term of 
the asymptotic expansion. The residue takes a value such that this divergence is exactly cancelled 
by limj,^_i I\2^^'^^\s). Explicit computation of (15) has been carried out for a = 2 and we have 
perfect concordance with the data depicted in the Figure 5. In this case, the coefficients c„ vanish, 
such that £A£-^(^)(s) = 0. 



6 Jacobi fields, kink instability and resonances 

In the C^^s^ topological sectors, things are more difficult. As critical points of the energy functional, 
the topological defects must satisfy the Euler- Lagrange equations 

^ = 201 (40? + 2a(l + a)^l - l) ^ = a02 (4(a + 1)0? + 2a0^ - l) . (16) 

\i a — 2 and cr = |, we are able to obtain all the kinks in these sectors of the model because the 
associated dynamical problem is integrable. 

• For a — I'm. the C^^^^ topological sectors, we find 

1 _ (-1)° (-1)^ sinh2(x + a) ^ 

<^TK2\x- h\ - 2 cosh2(x + a) + 6'^ + 2 cosh 2(a: + a) + " 

• For (T = |, in the C^^-j topological sectors we have the family of solutions: 

4>TK2\x-h\ = (-1)" sinh&sinh(x + a) ^ 

' cosh^(x + a) + 2 cosh h cosh (a; + a) + 1 



^ (^_-^^^ - sinh(3; + a) 

cosh^ (x + a) + 2 cosh h cosh (x + a) + 1 



62 
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For generic cr, only the TKl' kinks are known explicitly in C^^^: 

^TKi'{x) = (-1)"— ^tanhW ^{x + a) 
yza V ^ 

Defining z = -s/^x, we write the Hessian for the TKl' kinks as 

^"2!^ ' 0" _^ + 4-6sech2^ ; ■ ^ 

i7 is a diagonal matrix differential operator that includes as the J722 matrix element the very well 
known Posch-Teller Schrodinger operator, with Cq = 4 and uq — 6. Jn is also of Posch- Teller type 
- cl — A, uo — and deserves a closer analysis. The eigenvalues ujn — 4: — [k, — {n + |)]^, 

1^ = \J '^''^t'^^ + |"; ^ = 0, 1, 2, • • • , I[k — |] — 1 belong to the discrete spectrum of Ju. Thus, 

Spec J = {4 - - (n + U {kl + 4} U {0} U {3} U {kl + 4}^,^^ ■ 

The phase shifts for the continuous spectrum are 

3A; 

5{k]) = S+{ki) + 6-{ki) ; S(k2) — — 2arctan- — 

2 — k2 

The reflection coefficient associated to the Schrodinger operator Jn is not zero if cr 7^ j(j4^i)_4 
with j a positive integer number greater than or equal to 2. Therefore, there are also in general 
even and odd phase shifts, 

5^f/c,)- Cretan (^ ^"^(^(^^^^-^(^^^^ 

to be read from the transmission and reflection 

^ T{l + K-tk^)T{l-K-ik,) ^ T{^K-tk,)T{l-K^-tk,)T{tk,) 

^ r{l-iki)T{-zki) ' ^ T{l + K)T{l-K^)r{-iki) 

coefficients. 

The physical and geometrical information encoded in the spectrum of is as follows: 

• There are two bound states in the tangent direction to the TKl' kink in C^%y The first 
of them is the translational zero mode and the second one is interpreted as usual: a scalar 
boson "polarized" in internal space in the 6*2 direction is captured by the kink. The scattering 
eigenstates correspond to ci-polarized scalar bosons moving in the kink background. 

• There are I[k — ^] — 1 bound states in the orthogonal direction to the TKl' kink in C'^y 
The associated eigenvalues are positive if cr > 2 and thus the TKl' kink is stable in this 
regime. 
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When cr < 2, there are some negative eigenvalues in Spec(j7') and the TKl' kink is unstable 
against the associated fluctuations. The number of negative eigenvalues and the degree of 
instability increase with decreasing a: if cr e ( ^.^^._^^^_^ , the spectrum of Ju contains 

J — 2 negative eigenvalues and there exist j — 2 instability directions around the TKl' kink 
in the C'^^) sectors. For an analysis of kink stability in multi-component scalar field theory, 
see [33]. 

At the critical values a — there is a jump in the number of negative eigenvalues and 

the Jacobi field 



V'j(g) = (l-g')2Fi[i + Jj(j + l) + i|-\AXj + l) + i,3,i(l + g)] ; g = tanh^ 



becomes a bound state belonging to the discrete spectrum of Ju. The Morse index theorem, 
[33], tells us that there is a continuous family of TK2' kinks in the C^.^.) sectors if o" = j^j^^yri- 
The members of this family are the flow lines of the gradient of a second superpotential W 

that has two critical points at 0y± = ± ^ — 62 and j — 2 branching points: the zeroes 
of '0j(x)e2, see [33] . 




Figure 6: Graphics of il)j{q) for j = 2,j = 3,j = 4,j = 5 (left), for non-critical values of a close to the 
critical ones (right). 

Therefore, the TKl' kink gives rise to a bona flde eigenstate of the Hamiltonian only if cr > 2. 
Here we shall not give AM(TKl'), because the computation is identical to that developed for 
AM(TKl) if cr > 2 but gives rise to a complex mass if cr^ < 2; the TKl' kink state should be 
interpreted as a resonance in this regime. 



Appendix: mode number cutoff regularization method 

In this Appendix we present a derivation of the basic equation (6) in Section §3. Given the /C and 
V Schrodinger operators 
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respectively ruling the small fluctuations around the vacuum and the kink, the renormalized zero 
point energy is: 



Here, = ^:^|^/)v is a constant, and v'^ — V{x) = ^^I^ir gives rise to a potential well V{x), which 
rapidly decreases to its asymptotic value: lim V{x) = 0; note that, as in the main text, we 



use dimensionless eigenvalues and recover the dimension in the flnal response through the overall 
constant of (18). The spectrum of both fC and V is continuous and some regularization is needed to 
compute the infinite sums in the Casimir energy (18). Over the last few years it became clear that 
the correct regularization procedure is based on a cut-off in the number of fluctuation modes, a 
number that must be taken to be "equal" in the vacuum and kink sectors, see [12]. More recently, 
controversy arose about whether an energy cut-off regularization can yield the correct answer, see 
[34] and [35]. This is a subtle matter indeed, but all the evidence leads us to rely on the mode 
number regularization method. 

Before deriving (6), we outline the set-up pointed out in [14] that underlies our approach below. 
The problem is discretised by confining the system to a box of length 2L, which is much greater 
than any length scale characteristic of the potential term or of any of its bound states, and the 
limit L — > oo is taken at the end of the analysis. In order to specify the boundary conditions we 
distinguish between (a) eigenstates with eigenvalue greater than v^, giving rise to the continuous 
part of the spectrum, and (b) those with eigenvalue lower than or equal to v"^, corresponding 
respectively to bound or half-bound states. Fictitious boundary conditions ip{±L) = are imposed 
on the eigenstates of type (a), but we can safely anticipate the limit L — > cxo for states of type 
(b) and use for them the eigenfunctions of the actual bound or half-bound states respectively 
decrease exponentially or go to a finite constant at long distances. In any case, we stress the fact 
that the effects induced by the boundary conditions disappear when the L — > cx) limit is taken. 
Nevertheless, we shall approach first the problem by focusing on operators /C without half-bound 
states. This particular case (tackling operators which involve potentials without refiection, such 
as those in the kink and soliton of the well known (p'^ and Sine-Gordon model) will be dealt with 
later. 

Apart from restricting the spectra of V and /C 



to those eigenfunctions of type (a) satisfying Dirichlct boundary conditions i/j{±L) = 0, we shall 
only consider symmetric potentials V{x) = V{—x) which allows us the choice of odd, ip{x) = 
—ijj{—x), and, even, ijj{x) = iJj{—x), eigenstates complying with the boundary conditions. This 
breaks the degeneracy that arises by imposing periodic conditions on eigenfunctions with non 
defined-parity (plane waves) , and hence skips some problems arising when the number of bound 
states is odd. 

A brief description of the eigenvalues and eigenfunctions of V and /C in a finite box follows. 
1. Spectrum of V 

• Half-bound state. 




(18) 



Vip[k;x] = cu'^ {k)il)[k; x] 



K:ij[q;x] = uj^{q)ip[q;x] 
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The constant function, 

ipo{x) = constant , 
is an eigenfunction of V with eigenvalue f ^ that do not grow at a; = ±00. 

• Odd eigenfunctions 

If n~ is a natural number strictly positive and the wave vectors k^- satisfy 

kn- ■ L — Trn~ , n" e N"*" , 

then the odd functions, '^„-(a;) = —ipn-^—x), 

i>n- {x) = An- ■ sin(/c„- • x) 

belongtothespectrumofV with eigenvalue = A;^_+v^. Since {kn-+i — kn-)L — 

TT the density of odd eigenstates is: p~{k) — ^. 

• Even eigenfunctions 

If n'^ is a natural number and the wave vectors kn+ satisfy 

kn+-L^ 7r(n+ + i) , n+ e N 
then the even functions, ipn+ix) — ipn+{~x), 

1pn+ix) = An+ ■ COs{kn+ ■ x) 

belong to the spectrum of V with eigenvalue a;^(A;„+) = k'^+ + . Now, p^ik) — ^■ 

2. Spectrum of /C 

• Bound states: 

We assume that there exist / bound states with uP' < v"^. In particular we denote 
by l'^ the number of symmetric bound states and by l~ the number of antisymmetric 
states, / = + l~ . If Lof = v^, the highest bound state becomes a half-bound one and 
contributes with a weight of s; = ^ to the mass quantum correction. 

• Odd eigenfunctions. 

There are odd eigenfunctions which asymptotically are of the form: 

-J. 

V^„-(,t) ~ At- ■ sin[g„- • a; + 

The odd phase shifts 5' are defined in terms of the transmission and reflection 
coefficients, see Section §3, but now the boundary conditions require: 

•L + r(g„-) =7rn- , e N . (19) 

Unitarity, transparency at g = 00, and, continuity of the wave function at g = 
(threshold) allow us to set the phase ambiguity as in [14] : 

lim 5~{q) — , lim 5~{q) — 7rl~ 

q-*oo g— »0+ 
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The question arises: what is the minimum value n = tIq? For very large L, q^- is 
very small and (19) becomes: 



;L + (5'-(0+)) + r(0+) =7rno 



(20) 



Here, 5'~(0+) is the hmit of the left-hand derivative of the antisymmetric phase shift 
when q goes to zero from the right. The minimum value of n~ — Uq is therefore 
— l~ . Since {qn-+i — Qn-)L + 5~(g„-+i) — S{qn-) — tt the density of antisymmetric 
states is: 

, , L 1 dS-(k) 

TT TT ok 

Even eigenf unctions. 

There are even eigenfunctions which asymptotically are of the form: 

i/jn+ {x) ""ci^ An+ ■ cos[g„+ • X + 6^ {qn+)] 

The phase shifts S'^{qn+) are defined in terms of the transmission and reflection coeffi- 
cients, see Section §3, and the boundary conditions require: 



qn+ ■ L + 5+(g„+) = 7r(n+ + i) 



n+ e N 



(21) 



Again we set the phase ambiguity by invoking Levinson's theorem as in [14]: 



lim 5+{q) = 

q— >oo 



lim S+{q) = 7r(Z+ - |) 



where is the number of even bound states. 
Prom 

g„.-(L + 5'+(0+)) + 5+(0+)=7r« + i) , 
we read the minimum possible value of q^+: Uq — l^. Now, 

TT TT ok 

We stress that due to the parity of the eigenfunctions there is no need to consider 
negative values of n^. 

The regularized version of formula (18) reads: 



K 



i=l 



+ 



+ 



hm 

~Y 

hm 



N 



No 



.n+=l+ 
N' 



n+=0 



n =1 



n-=0 



+ 
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The first row collects the contributions of the I = + bound states of A^. In the other two 
rows, the contributions of even and odd eigenfunctions in the continuous spectra are accounted 
for. Note that we have added the contribution of the constant mode of V to the odd sector because 
it corresponds to the lower eigenvalue of V. We stick to mode number cut-off regularization, i.e. 
we impose N = Nq and N' = N^. This choice balances the number of states of /C and V involved 
in the above formula. We have distinguished the number of symmetric and antisymmetric states, 
N and N', although we shall consider either N' ^ N + 1 or N' ^ N. 

To measure the contribution of /C with respect to the contribution of V to the Casimir energy 
mode by mode, we write the second row in the form: 

N / X Z+-1 



hm 

~Y 

hm 



.n+=0 
N 

-E 

. n+=0 



n+=0 



1+-1 



n+=0 



L — *oo 



hm 



Jo dk 



A similar process for the odd eigenfunctions gives: 



hm 
~Y 

hm 



N' 



L — '■oo 



hm 



E 

.n-=0 

N' 

-E 

. 71- =0 

oo 



q2_ _^ ^2 



V n 



kl. + 



E V^n-+-^ 



n-=0 



-E v^n-+^'+^(^; 



n-=0 



dk dL}{k) 



TT 



5-{k)^l-v 



The sum of all contributions plus a partial integration. 



hm 



I 



hm 



2 
hm 



-l+v 



^s^ikMk) 



I'v 



hm 
2^ 



5-{k)uj{k) 





oo 



hm 
27 
hm 
2^ 



dk 



dk 
d6-{k) 
dk 



u{k) 



and the asymptotic behaviour of the total phase shift. 



fc— >oo 1 

5^{k) - - 



dxV{x) co^{kx) 



fc— >oo 1 



dxV{x) s\T?{kx) 



show that: 



hm sr^ hm , , hm ^ , , , hm , hm ^ , , 



hm 
' 2^ 



/•oo 

/ dxV[x)\ 
Jo 



27r ' ' 2 
cos^(A;x) + sin^(A;x)] + 



2n 



hm r 

Jo 



dk 
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From the ID Levinson theorem we notice that the two last terms in the first hne of the above 
formula cancel each other but the two preceding ones leave a contribution: — |. Therefore, the 
formula for the kink Casimir energy reads: 



2 



E 



V 1 



dk 



27r 



(22) 



where {V{x)) = V{x)dx — 2 V{x)dx is the expectation value of the potential term and the 
contribution of the continuous spectrum is encoded in the total phase shift: 6{k) = S'^{k) + S~{k). 
We call attention to the contribution — ^hbh which tells us that the contribution of the mode with 



k 



- wave vector in SpecV is half the contribution of a bona fide bound state -hence, the name-. 
The generahzation of this approach to tackle the presence of a bound state with k — and 
= v"^ in Spec/C (half-bound state) is direct. The Levinson theorem. 



5_(0+) = /_7r 



5+(0+)=Z+7r 



TT 

2 



(23) 



works as above, but there is now a half-bound state that contributes with | to both the number 
/ = and the mass quantum correction Ai^^, see [13]. Thus, we write (23) as 5(0+) = 

Ztt — I = n^TT, and, bearing in mind the previous observation, the above calculations reproduce 
the results obtained in Reference [12] . The tricky point to realize is how the contribution of the 
threshold state turns (23) into 5(g+) = n(,7r, with rib the number of bound states with o;^ < v^. 
In the kink mass formula, the contribution of the half-state would cancel out the term — |, and 
everything is okay. For this reason, the formula shown in [12] is valid only for operators involving 
potentials without reflection. A final expression involving both of these cases provides the first 
row of formula (6) 



hm 



i-i 



^ V 1 

i=l 



dk 



27r 



(24) 



where si ■ 
non-zero. 



i if we deal with a reflectionless potential but = 1 if the reflection coefficient is 
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